Dynamical process for switching between the metastable ordered magnetic state and 
the nonmagnetic ground state in photoinduced phase transition 
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We propose a dynamical mechanism of the two-way switch- 
ing between the metastable state and the stable state, which 
has been found in experiments of photoinduced reversible 
magnetization and photoinduced structural phase transition. 
We find that the two-way switching with a non-symmetry 
breaking perturbation such as illumination is possible only in 
systems with appropriate parameters. We make it clear that 
the existence of two time scales in the dynamical process is 
important for the two-way switching. 



Dynamics of the switching between the metastable 
state and the stable state is one of the most interesting 
topics of phase transitions. Photoinduced phase transi- 
tion (PIPT) has currently attracted interests in the con- 
densed matter physics Min- In spin-crossover complex 
several such phenomena have been reported [pImIH] • PIPT 
is accompanied by a structural phase transition as ob- 
served in TTF-CA H or by a magnetic phase transition 
as observed in Co-Fe complex plpf. Because the energy 
level of the low spin state of such a complex gets cross to 
that of the high spin state, the spin state can be changed 
between these states by stimulus due to changes of tem- 
perature, pressure, photoexcitation, etc. Among them, 
the cobalt-iron prussian blue analogs [g]||] have been re- 
ported to be able to switch from the paramagnetic state 
to the ferrimagnetic state (or vice versa) by visible (or 
near-IR) illumination. 

In order to study the switching mechanism, we consider 
a system which has a long-lived metastable ferromagnetic 
excited state. The transition process from the metastable 
magnetic state to the stable paramagnetic state can be 
realized by the nucleation process because the bulk free 
energy of the stable state is lower than that of the initial 
metastable state || . Hereafter we refer to the metastable 
state as MS and to the stable state as SS. On the other 
hand, switching from SS to MS conflicts with the ther- 
modynamic stability and it is considered to be difficult 
to realize as far as we apply only symmetric disturbance, 
e.g., illumination, although it is easy by using a uniform 
field which causes the symmetry breaking. Our interests 
is how the two-way switching is realized in a symmetric 
applied field. In our previous study E3], photon's effect 
is attributed to the change of renormalized parameters 
(D,T), where we chose these parameters in the phase 
diagram from a viewpoint of phenomenology. 



In this letter we consider the condition for the two-way 
switching in terms of the transition probabilities among 
microscopic states. We find that the two-way switch- 
ing cannot be realized when we introduce only a spin- 
pumping effect by photon. In addition to the pumping 
effect of photon, a suitable relation between transition 
probabilities among elementary processes is necessary to 
realize the two-way switching. Here we will extend the 
Glauber dynamics |l3| and introduce a more complicated 
relaxation process where the switching (SS — » MS) can be 
realized. Using this model we will study characteristics 
of the system which can show the two-way photoinduced 
transition. 

As a typical model which shows the first-order phase 
transition, we adopt the Blume-Capel (BC) model JTI 12 
in the simple cubic lattice, 
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where Si=H or and (i,j) denotes the nearest-neighbor 
pairs. Here Si = ±1 represents the high spin state with 
the magnetization ±1, respectively, and Si — represents 
the low spin state which is nonmagnetic. The excitation 
energy from the low spin state to the high spin state is 
D(> 0). We take J(> 0) as a unit of energy and set 

*fe = l. 

Here we adopt a stochastic dynamics with multi-time 
scales as follows. Let P(Si ■ ■ ■ S k ■ ■ ■ SjVj t) be the proba- 
bility of the state {Si ■ ■ ■ Sk- ■ ■ Sn} at time t. Using the 
transition probability per a unit time Wk(Sk — * S' k ) at 
the k-th site (Fig. [j] (a)), the master equation is generally 
given by 

±P(S 1 ---S k ---S N ,t)= (2) 

at 

-£ £ W k (S k ^S' k )P(S 1 ---S k ---S N ,t) 
+ J2J2 W k (S' k ^S k )P(Si---S k ---S N ,t). 

k S k =£S' k 

Here, we assume that the time scale of the transition 
between S = ±1 is much faster than that between 
5 = and \S\ — 1, because the former is a simple spin 
flip process while the latter is accompanied by a struc- 
tural change. Representing this situation, we introduce 
two transition processes: The first one is the standard 
Glauber type evolution, 
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where S fc = ±1, 0(Fig. | (a)) and y = /3J£ Nr ,S;. 
Here X)nn denotes the summation over the nearest neigh- 
bor sites of k. We express this process by 4zP(t) — 
LG(t)P(i). The second one is a process of spin flip, Lq, 
(Fig. 0(b)), 
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We study a combination of these two processes, 
L(t)= PQ L G (t)+psLs(t). 



(5) 



Here the parameters pc and Ps (pg + Ps = 1) describe 
the relative time scale of the processes. When pg/ps is 
small, the relaxation among the high spin states (5 = 
±1) occurs more quickly than that between high and low 
spin states (\S\ = 1 and 0). Because both Lq and Lg 
satisfy the detailed balance, it is warranted that with 
any combination of pg(^ 0) and ps any initial state goes 
to the equilibrium state. Here we think that the ratio of 
PG and ps is given as an inherent property of individual 
materials. Under the process Ls, the pumped spins tend 
to align before the system relaxes back to the state of 
5 = 0. When the adequate balance of these two processes 
is realized, we expect that the reversible switching of the 
magnetization by illumination is possible. 
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FIG. 1. Transitions among states. The circle is the fcth 
site. | and J. in the circle denote up spin state and down spin 
state, respectively. Blank circle denotes a nonmagnetic state, 
(a) transition process among all the possible local states, (b) 
transition process between the high spin states S — ±1. (c) 
transition process by pumping. A in the circle denotes any 
state (5fc = l, —1 or 0). 



Next we consider a process Lp for illumination. During 
the illumination, the pumping process is included. Pho- 
ton excites a site which is either nonmagnetic or magnetic 
to a magnetic site. It is a natural assumption that the 
transition to up spin or to down spin has an equal prob- 
ability and we take the transition probabilities (Fig. |l| 
(c)), 



W k {S k - ±1) = -. (6) 

The time evolution during the illumination is given by 



L(t) = p G L G (t) +p S L S (t) +p P L P (t), 



(7) 



where pg + Ps + Pp = 1< The magnitude of pp repre- 
sents the strength of the excitation. Indeed in the ex- 
periments [|2|,|3| , two kinds of illuminations with different 
frequencies are used to realize the switching of magneti- 
zation. 

We look for a suitable set of po, Ps and pp with which 
the system can be switched between MS and SS. Let us 
study the dynamics of the system first within a mean 
field theory (MFT). We will then check the stability of 
the results of MFT by a Monte Carlo (MC) method. 

In MFT, Eq. (||) is reduced to an equation of pi(t), 
p_i(i), and po(t), which are respectively the probabili- 
ties for up spin state, down spin state, and nonmagnetic 
state on a site at time t. Furthermore, the equations of 
pi(i) and p_i(i) are reduced to an equation of the mag- 
netization, and the equation is found to be expressed in 
the form of the Van Hove's phenomenological equation. 

d . , N , , , df(m, t) , . 

- m =p 1 (t)-p_ 1 (t) = -b^-J-L. (8) 

at am 

This f(m, t) is regarded as the effective free energy, which 
is useful to visualize the thermodynamic force at a given 
time t. We obtain f(m,t) by integrating the right hand 
side of Eq. (g) as to m. 

f(m,t) = -Jzm 2 ~ Pg— ln[2cosh(/3Jzm) +exp(/3D)] 

-p s -{ ln[cosh(/3 Jzm)] (1 - p (t = 0) exp(-(p G +pp)t)) 



PG 



-(1 — exp(pc +pp)i)(ln[cosh(/3 Jzm)} 



Pg +Pp 

ln[2 cosh(/3 Jzm) + exp(/?£>)])}. 



(9) 



In the case of no illumination (pp = 0), the above free 
energy gives the equilibrium free energy / oq (m) in the 
limit t — > oo regardless of the parameters pg(t^ 0) and 
ps- In Fig. | / cq (m) for (D,T) = (3.2,0.6) is shown. 

Because of the existence of metastable valleys, some 
initial states may relax to the MS instead of SS. In the 
case of ps =0 (pc = 1), f(in,i) is independent of time 
and takes the form of / eq (m) in Eq.ffl). When the value 
of magnetization is larger than the critical value m c (= 



±0.5557), the state relaxes to MS. Thus, in order to make 
the system to relax to MS after the illumination, the 
magnetization in the illumination must increase beyond 
m c . It is found however that the magnetization of the 
stationary state in the illumination is much smaller than 
m c for any value of pp. Therefore it is impossible to cause 
the transition SS — > MS when ps = 0. 

Next we consider the case of ps ^ 0, where the process 
of spin ordering is expected to be enhanced. In spite of 
the case of ps = 0, the form of f(m, t) depends on t and 
Po (t = 0) . In Fig. pi we show a kind of phase diagram 
of the initial values (po, to) for the cases of (a)pc/Ps = 
0.667, (b)p G /p s = 0.163 and (c) p G /ps = 0.429. Because 
of the condition pi+p_i+po = 1 and < pi,P—i,Po < 1, 
Po + m < 1 must be satisfied. The open circles denote the 
boundary between the region of (po, to) from where the 
state goes to MS and the region from where the state goes 
to SS. There the values (po, to) of the stationary state 
of the illumination process (pp ^ 0) are also plotted by 
closed circles for various values of pp keeping the ratio 
of pg/ps- That is, in the case of pp = 0.5 and pg/ps — 
0.667, p G = 0.2 and p s = 0.3. 

If a closed circle for a value of pp is located below the 
critical line linking open circles, the initial state with the 
values (po, m) denoted by this closed circle relaxes to 
SS (para). By the illumination the system becomes to 
the state denoted by this closed circle, and then after 
the illumination it relaxes to SS. In the case of (a) the 
point (po, to) of the stationary state with any value of 
pp is located below the critical line as shown in Fig. 3 
(a), and thus all states after the illumination give the 
initial state for SS. On the other hand, if a closed circle 
for a value pp is located above the critical line, the state 
relaxes to MS. By this illumination the system becomes 
to the state denoted by this closed circle, and then after 
the illumination it relaxes to MS (ferro). In the case of 
(b) the point of the stationary state with any value of pp 
is located above the critical line as shown in Fig. 3 (b), 
and thus all states after the illumination give the initial 
state for MS. In order to realize the two-way switching 
(MS=^SS) the system has to have both cases for different 
values of pp. That is, the closed circles are located on 
both sides of the critical line as shown in Fig. 3(c). 
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FIG. 3. Critical line linking the open circles in the relax- 
ation process (pp = 0) and the values (po, to) of the stationary 
state in the illumination are also plotted by the closed circle 
for several values of pp. (a)pG/ps = 0.667 (d)pg/ps = 0.163 
(c)pg/ps = 0.429 

With the present observation we find that the two-way 
switching is not possible in every system but in only lim- 
ited materials which have a peculiar property that the 
boundary (open circles) crosses with the line of closed 
circles. On the other hand, we also find that such situa- 
tion can be actually exist in some appropriate condition. 

Now we demonstrate how the switching of the magne- 
tization is realized. According to the results in Fig. 3, we 
choose the ratio pg/ps =0.429. As an initial state we take 
the equilibrium state of the parameter (D, T) = (3.2, 0.6) 
where to — and po = 0.99. We make the system in the 
illumination of pp = 0.7 for a period to < t < t\ (pro- 
cess I), and then we turn off the illumination and let 
the system relax, i.e. pp = for a period t\ < t < tn 
(process II). The change of the effective free energy 
f(m,t) in the process I is shown in Fig. (a). At the 
end of the process /(to, t) almost reaches the form of 
its stationary state where /(to, t) has two minima and 
to = ±0.18. In the second process f(m,t) changes as 
shown in Fig. || (b). Because the final state of the pro- 
cess I (po, to) = (0.11, 0.18) locates in the region to go to 
MS (ferro), the system relaxes to the metastable state. 
Thus we have the switching SS — > MS. We see that the 
magnetization increases with time in the process II. 



FIG. 2. Mean field free energy at the equilibrium state for 
(AT) = (3.2,0.6). 
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FIG. 4. Time evolution of f(m, t) fixing the ratio 

Pg/ps = 0.429. (a) process I, to < t < ti (pp — 0.7), (b) pro- 
cess II, ti < t < tn, (c) process III, tn < t < ini (pp = 0.9), 
and (d) process IV, tin < t < irv- The thick arrows show the 
change of the magnetization. The change of f(m, t) is shown 
for t a < t < i a + 20 (a=0, I, II and III) which corresponds 
to the first 1/5 of the process because f(m, i) changes little 
after t a + 20. Here the sign of magnetization is not essential. 

Next, let us consider a procedure for the switching MS 
— > SS. Now we illuminate the system with another light 
of pp = 0.9 for a period tn <t < tin (process III), where 
f(m,t) changes as shown in Fig. H (c). The f(m,t) 
at the final state (t = tin) is parabolic and the system 
has almost no magnetization. Thus after turning off the 
light the system relaxes to the stable paramagnetic state 
as shown in Fig. H (d) (process IV). The time evolution 
of magnetization of the whole processes is shown in Fig. 
^ (a). In the process III we could choose pp to be 0.4 
instead of 0.9 to realize MS— >SS because the closed circle 
for pp = 0.4 is also located below the critical line. 

We have shown that the two-way switching is real- 
ized by the processes I ~ IV in MFT. In order to check 
whether the switching process is stable against fluctu- 
ation, we study the processes by a Monte Carlo (MC) 
method. Of course quantitatively MC and MFT gives 
different results. But, qualitatively the same features of 
the dynamics as obtained in MFT are also found in MC. 
Namely, when ps = 0, the switching SS — > MS is not 
observed for any value of pp. On the other hand we 
can realize the two-way switching by tuning the param- 
eters as pg/ps = 0.289, pp = 0.096 for the process I 
and pp = 0.24 for the process III. The time evolution of 
magnetization per site is shown in Fig. pj (b) ( simple 
cubic lattice with 1000 sites in periodic boundary condi- 
tion). Strictly speaking, in MC the true stationary state 
of the process II is the paramagnetic state, but the relax- 



ation MS — ► SS is very slow and thus practically the same 
behavior as in MFT is reproduced. We confirmed that 
this switching of magnetization can be repeated reliably, 
and we conclude that the two-way switching provided by 
the processes I - IV is stable even in systems with short 
range fluctuations. In the MC simulation, the switching 
dynamics depends very sensitively on the ratio of param- 
eters as well as in the MFT case. 



(a) 



process I 

/ 



I ;\ process IV 

, process 111, 



200 300 400 

ti lime t[[ tin 



"(b) 



0.8 - 
0.6 - 



proce ss II 



process I 



process III 



it 



process IV 



1 10 s 


2 10 5 


3 10 s 


4ltf 


5 Ltf 


to 


ti 


MCS tn 


tin 


tiv 



FIG. 5. (a)Time evolution of magnetization per site for 
the parameter set pg/ps = 0.429. pp = 0.7 and pp = 0.9 are 
adopted as the first and the second illumination respectively 
in the mean field theory, (b) Time evolution of magnetization 
per site for the parameter set pg/ps = 0.289. pp = 0.096 and 
pp = 0.24 are adopted as the first and the second illumination 
respectively in the Monte Carlo method. 

In the evolution Lp, we considered that illumination 
causes the low spin state to be excited to the high spin 
state. However, illumination may also cause the high 
spin state to change to the low spin state, which disturbs 
the development of magnetization and gives a negative 
effect on the switching SS — > MS. As the extreme case 
we consider that illumination brings one of three states 
with the equal probability regardless of the initial state, 
i.e., W k (S k ->• S' k ) = | for any S k and S' k . Although it 
is difficult to grow a magnetic droplet in this case, we 
found that the switching (SS ^ MS) is still possible. 

We studied dynamical processes for photoinduced two- 
way switching between the ordered state and the disor- 
dered state in a master equation approach. We have 
found that a peculiar condition is necessary between the 
transition probabilities. That is, only when the transi- 
tion probability between the high and low spin states and 
that between the high spin states satisfy a special ratio, 
the two-way switching by using different illuminations 



is possible. This implies that only materials which sat- 
isfy such conditions can show the photoinduced switching 
phenomena. 
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